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Determination of the Upper Critical Field of a Single Crystal LiFeAs: The Magnetic 

Torque Study up to 35 Tesla 
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We report on the upper critical field B C 2 of a superconducting LiFeAs single crystal with 
T c ~1 6 K, determined from magnetic torque measurements in dc-magnetic fields up to 35 T and 
at temperatures down to 0.3 K. B C 2 at 0.3 K is obtained to be 26.4 T and 15.5 T for the applied 
field fl a || ab and B a ||c, respectively. The anisotropy parameter F = B^2 /B^2 is ~ 3 at T c and 
decreases to 1.7 as T — > 0, showing rather isotropic superconductivity. While B C 2 is orbitally- 
limited for _B a \\c, the spin-paramagnetic effect is evident in the temperature dependence of 
B C 2 for B a || ab. 
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Since the discovery of superconductivity in 
LaFeAs(0,F) with T C = 26K, 1 ) a variety of related 
compounds containing FeAs-layers has been found to 
exhibit superconductivity. 2 - 1 The parent compounds 
iJFeAsO (i? = rare earth, "1111" system) and AFe 2 As 2 
(A = alkaline earth or Eu, "122" system) with the 
ZrCuSiAs- and ThCr 2 Si2-type structures, respectively, 
undergo antifcrromagnetic and structural transitions. 
The transitions can be suppressed by several kinds of 
doping effects 1 ' 3-6 - 1 or application of pressure, 7 ~ 9 ) and T c 
reaches ^56K in some compounds. 10 ~ 12 ) The magnetic 
long range order usually competes with superconduc- 
tivity, but the fluctuation likely plays a crucial role in 
the pairing mechanism of the Fe-based high-T c systems. 
It is also of interest that, around the optimal condition 
where T c shows its maximum, deviation from conven- 
tional Fermi-liquid behavior has been observed such as 
p^T, 13,14 ) anomalous Hall angle, 13 ) an enhancement 
of effective masses 15 ^ etc. These superconducting and 
normal-state features bear resemblance with those 
widely reported for strongly-correlated electron systems 
including cuprates and heavy fermion compounds. 

The title compound LiFeAs, categorized into the "111" 
system with CeFeSi-type structure, has distinctive char- 
acteristics: (i) the stoichiometric superconductivity with 
T c as high as ~ 17 K, 16 ) (ii) no experimental evidence for 
the magnetic/structural transitions, 17 ' 18 - 1 and (hi) single 
crystals with high quality (residual resistivity ratio up 
to 50). 19 ' 20 * 1 Therefore, LiFeAs provides a unique oppor- 
tunity to probe the intrinsic properties of the Fe-based 
high-T c superconductivity. 

Up to now, there are few reports on systematic mea- 
surements of the upper critical field B c2 of Fe-based su- 
perconductors to address the issue of the pair-breaking 
mechanism. This is mainly due to the fact that high- 
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T c superconductors including Fe-based systems generally 
have extremely high £? c2 . In many cases, accordingly, the 
low-temperature behavior is extrapolated from the high 
temperature data around T c , which may lead to mislead- 
ing conclusions. A precise determination of f? c2 over a 
whole temperature range could provide important clues 
to the pair-breaking mechanism of high-T c superconduc- 
tivity. 

Here, we present the first report on the whole temper- 
ature dependence of B c2 and its angular variation of a 
LiFeAs single crystal, using a magnetic torque technique 
with a 35 T dc resistive magnet. 

Single crystals LiFeAs were prepared from high-purity 
constituent elements in a ratio of Li:Fe:As = 2:1:2 by a 
self-flux method. 20 ' 21 ) The residual resistivity ratio is as 
high as 45, showing the high quality of the crystals. Tem- 
perature dependence of the dc magnetic susceptibility of 
LiFeAs was measured with an applied field B a of 1 mT 
for B a || ab, in a Magnetic Property Measurement Sys- 
tem (MPMS: Quantum Design) around T c . As displayed 
in the main panel of Fig. 1, the bulk superconductiv- 
ity was confirmed from the clear diamagnetic transition 
below T c ~ 16 K with approximately 16 % (97 %) of the 
superconducting Meissner (shielding) volume fraction for 
(zero-)field-cooled process. A tiny piece taken from the 
same crystal was mounted on a piezoresistive microcan- 
tilever with a small amount of grease (see, inset photo 
of Fig. 1), attached to a rotator probe, and inserted into 
a 3 He cryostat (Heliox: Oxford). The angle is the an- 
gle between the applied field B a and the c axis: 6 = 0° 
for B a || c and 9 = 90 for B a \\ ab. The magnetic torque 
was measured with a water-cooled dc resistive magnet 
in fields up to 35 T and at temperatures down to 0.3 K. 
Since the magnetic torque, proportional to M x B a , is 
very small at 6 = 0°, the torque is measured at 6 = 5° 
to improve the signal-noise ratio, which is specified as 
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Fig. 1. (Color online) Magnetization divided by applied magnetic 
field, M/B&, vs T of a single crystal LiFeAs around the transi- 
tion temperature, T c , in a field of B a = 1 mT applied parallel to 
the ab plane. Results measured in zero-field-cooled (ZFC) and 
field-cooled (FC) processes are displayed. A tiny piece cleaved 
from the same crystal was mounted on a microcantilever for the 
magnetic torque measurement as shown in the inset photograph. 
The field angle 9 is the angle between the c axis and the applied 
magnetic field B a . 
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B a || c. We confirmed that the difference of B c2 between 
the two angles is negligible from the angular dependence 
oiB c2 . 

Figure 2 shows the field dependence of the torque sig- 
nals of LiFeAs up to 35 T at fixed temperatures down to 
0.3 K for (a) B a \\ ab and (b) B a \\ c. One can see highly 
hysteretic behavior between the field - up and - down 
sweeps, which is more significant at low temperatures. 22 ) 
The hysteresis of the torque response, the irreversible 
curve, appears in the superconducting mixed state when 
the pinning force is strong enough to trap the flux lines. 
In this study, we define the upper critical field B c2 as 
the field where the irreversibility disappears. 23 ) To accu- 
rately determine B c2 by minimizing a drift of the torque 
signals, we made the small-loop measurements around 
B c2 at each temperature, as shown by the dashed curves 
in (a) and (b) for T = 0.3K, in addition to the large-loop 
ones between B a = and above B c2 . The definition of B c2 
is illustrated in the inset of Fig. 2(a), where the vertical 
axis represents the difference of torque signals between 
field-up and -down sweeps of a small loop. Note that the 
shapes of the small loops near B c2 for B a || ab and B a \\ c 
are quite different. It could be possible that the pinning 
mechanism is different for the two field orientations, al- 
though the detail remains to be clarified. 

Figure 3 displays thus determined B c2 of LiFeAs 
as a function of T for £> a ||a6 and B a \\ c. The B c2 
curves are consistent with the results obtained by spe- 
cific heat (C), ac-susceptibility (x ac ), and 75 As nuclear 
magnetic resonance (NMR) measurements using sam- 
ples of the same batch , 24 ) and are qualitatively sim- 
ilar to other results. 25,26 ) As indicated by the dashed 
curves, the data can be well fitted using a Werthamer- 



Fig. 2. (Color online) Magnetic torque signals of LiFeAs as a 
function of B a at various temperatures down to 0.3 K for (a) 
B a || ab and (b) B a || c (inset: magnified view). The solid curves 
represent the data obtained in field-up and -down sweeps be- 
tween B a = and above B c 2 , where the arrows indicate the sweep 
directions. Small-loop measurements around B c 2, as shown by 
the dashed curves at T = 0.3K, were performed at each tem- 
perature. In this study, B c 2 was defined as the field where the 
difference between field - up and - down signals in a small loop be- 
comes zero within experimental error as illustrated in the inset 
of (a). 

Helfand-Hohenberg (WHH) formula containing the spin- 
paramagnetic and orbital pair-breaking effects. 27 -' The 
fits give the Maki parameter a = 2.30 and 0.75 for 
B a || ab and B a || c, respectively, where T C = 15.5K is 
fixed. 28 ) Because the torque signal, ocMx£? a , dimin- 
ishes at low fields, it is difficult to unambiguously de- 
termine the initial slope dB c2 /dT\T=T c - Instead, the 
values of dB c2 /dT\r=T c are estimated to be 4.43 and 
1.44 T/K for B a \\ ab and B a || c, respectively, from the 
relation a — —0.52dB c2 /dT\T=T c - These values are com- 
parable to those obtained from ac-x- 24 ) The orbital 
critical fields B* 2 at T = are estimated to be 47.2 
and 15.3 T for B a \\ab and B a || c, respectively, from 
a dirty limit formula S* 2 (0) = 0.69 T c diJ c2 /dT| T=Tc . 27 ) 
The Pauli-Clogston paramagnetic limit B po ~ 1.84 T c 29 ) 
is 28.5 T. The Ginzburg-Landau (GL) coherence length 
£ is obtained to be = 4.64nm and £ c = 1.50 nm, using 
B* 2 ab (0) = $o/27rU(0)&(0) and B* 2 C (0) = $ /2< Qb (0) 2 , 
where $ = 27r?i/2e = 2.07 x 10" 15 Tm 2 is the flux quan- 
tum. At T = T C , the anisotropy of the effective masses, 
m* 6 /m*, is about 0.11, using B* 2 ab / B*J = (m*/m* ab ) - 5 . 
These superconducting parameters are summarized in 
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Fig. 3. (Color online)B c 2 vs T of LiFeAs for _B a || ab and B a || c. 
The dashed curves indicate fits to the data based on the WHH 
theory (see text). For comparison, data deduced from C, Xac, and 
75 As-NMR measurements are shown. 24 ' The inset shows the T- 
dependence of the anisotropy parameter T = / ■ The data 
at T c (•) is obtained from the WHH fit. 
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Fig. 4. (Color online) (a) Torque signals of LiFeAs as a function 
of B a at 0.3 K for several angles 9, where 9 = and 90° corre- 
spond to Bn || c and B a \\ ab, respectively, (b) 0-dependence of 
B c 2 determined at T = 0.3K. The solid curve is a fit by a GL 
formula (see text). 



Table I. Superconducting parameters of LiFeAs and 
KFe 2 As 2 , 30) for comparison, obtained from the WHH 
fits; a: Maki parameter, A so : Spin-orbit scattering parameter, 
B* 2 '- Orbital critical field, B po : Paramagnetic critical field, £: 
GL coherence length. 
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the Table I, together with those of a stoichiometric su- 
perconductor KFe2As2, 30 ' for comparison. 

It is interesting to note that B* 2 C {0) < B po < B* 2 ab (0). 
This relation results in the orbitally limited B C 2 for 
£> a || c, as indicated by Ag = oo, and strongly spin- 
par amagnetically limited B C 2 for B a \\ ab, as indicated by 
that B^{0.3 K)<B* 2 afc (0), and effectively reduces the 
B C 2 anisotropy at low temperatures. A similar trend, 
namely the trend that the weaker orbital effect for B a \\ab 
is partly compensated for by the spin-paramagnetic ef- 
fect to yield a reduced anisotropy, is clearly seen in a 
stoichiometric superconductor KFe2As2 with low T c of 
2.8 K 30 ) (see, Table I), as well as in the high-T c systems 
including (Ba,K)Fe 2 As 2 , 31 ' 32) Ba(Fe,Co) 2 As 2 , 33 ' 34) and 
"H"-typeFe(Se,Te). 35 - 38 ) 

Figure 4 shows (a) the torque signals vs B a at 
0.3 K for several angles, 9, and (b) B c2 vs 9 at 



0.3 K. The solid curve in Fig. 4(b) is a fit to 
the B C 2(9) data using a Ginzburg-Landau (GL) for- 
mula, 39 ) B c2 (0) = B c c2 / {cos 2 9 + Y- 2 sin 2 6») 5 . The fit 
yields B c c2 = 16.3 T and r = 1.63. Note that the GL for- 
mula is based on the orbital effect whereas the spin- 
paramagnetic effect is crucial for B a \\ab in the present 
case. The quality of the fit is therefore rather limited as 
indicated by the deviations of the data points from the 
fit curve. The inset of Fig. 3 shows the T-dependcncc 
of the anisotropy parameter T, defined as L = B^ 2 / B^ 2 . 
As temperature decreases, T decreases from T = 3.1 at 
T c (from the WHH fit). At 0.3 K, B^| = 26.4T and 
B^ 2 — 15.5 T, which give T = 1.7. The low-temperature 
variation of T is similar to those observed in the "122" 
and "11" systems , 31 ~ 38 ) but the value of T in LiFeAs is 
slightly larger. In the "122" and "11" systems, T = 2~3 
at T c and T appears to approach ~ 1 at K . 31 ~ 38 ) This 
is markedly different from the results reported for the 
"1111" system. In NdFeAs(0,F), for example, T is 9.2 
at T c 40 ) (no report for V as T —> because of the large 
B c2 ). The larger value of T in the "1111" system can 
be ascribed to the more two-dimensional Fermi surface 
structures than those in the "111", "122" and "11" sys- 
tems, as expected from band structure calculations . 41 ) 
Similarly, the slightly larger T of LiFeAs than those of 
the "122" and "11" systems might be due to the more 
two-dimensionality of LiFeAs . 41 - ) In KFe2As2, however, 
there is an exceptionally large anisotropy (r = 6.8 at T c , 
see Table I) for the "122" system , 30 ) which is, on the 
other hand, in accord with the large resistivity anisotropy 
Pel Pab ~ 40 . 30 - ) In contrast, p c /p a b is reported to be as 
small as 3.3 in LiFeAs, and 2 - 5 in Ba(Fe,Co) 2 As 2 42) 
To conclude, we have performed high-field magnetic 
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torque measurements of a LiFeAs single crystal up to 
35 T, and determined the T- and ^-dependence of B c2 
down to 0.3 K. The anisotropy parameter T = / B^ 2 is 
slightly larger than typical values found in the "122" and 
"11" systems, but is still small, indicating quasi-isotropic 
superconductivity. Its temperature dependence is simi- 
lar to those observed in the "122" and "11" systems. 
The detailed analyses show that, while B c2 for B a \\ c is 
limited by the orbital effect, the weaker orbital effect for 
B a || ab due to the mass anisotropy is partly compensated 
for by the spin-paramagnetic effect leading to a reduced 
anisotropy at low temperatures. 
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